Observational constraints on spatial anisotropy of G from orbital 

motions 

Lorenzo lorio^ 

Ministero dell'Istruzione, deirUniversita e della Ricerca (M.I.U.R.)-Istruzione 
Fellow of the Royal Astronomical Society (F.R.A.S.) 
International Institute for Theoretical Physics and High Mathematics Einstein- Galilei 
Permanent address for correspondence: Viale Unita di Italia 68, 70125, Bari (BA), Italy 

lorenzo . iorioOlibero . it 



Received 



accepted 



- 2 - 



ABSTRACT 



A phenomenological anisotropic variation AG/G of the Newtonian gravita- 
tional coupling parameter G, if real, would affect the orbital dynamics of a two- 
body gravitationally bound system in a specific way. We analytically work out 
the long-term effects that such a putative modification of the usual Newtonian 
inverse-square law would induce on the trajectory of a test particle orbiting a 
central mass. Without making any a-priori simplifying assumptions concern- 
ing the orbital configuration of the test particle, it turns out that its osculating 
semi-major axis a, eccentricity e, pericenter and mean anomaly Ai undergo 
long-term temporal variations, while the inclination I and the node i? are left 
unaffected. Moreover, the radial and the transverse components of the position 
and the velocity vectors r and v of the test particle experience non-vanishing 
changes per orbit, contrary to the out-of-plane ones. Then, we compute our the- 
oretical predictions for some of the major bodies of the solar system by orienting 
the gradient of G{r) towards the Galactic Center and keeping it fixed over the 
characteristic timescales involved. By comparing our calculation to the latest 
observational determinations for the same bodies, we infer AG/G < 10^^^ over 
about 1 au. Finally, we consider also the Supermassive Black Hole hosted by the 
Galactic Center in Sgr A* and the main sequence star S2 orbiting it in about 16 
yr, obtaining just AG/G < 10~^ over 1 kau. 

Subject headings: Experimental studies of gravity; Experimental tests of gravitational 
theories; Modified theories of gravity; Orbital and rotational dynamics 
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PACS: 04.80.-y; 04.80.Cc; 04.50.Kd; 96.25.De 



1. Introduction 



The possibility that the Newtonian coupling parameter G may experience 
macroscopic spacetime variations ranging from laboratory to cosmological scales 



has been investigated both theoretically ( Sciamal Il953 



19621: iLong] Il974l : iGershtevn fc Gershtevn 



1988 



Linde 



Brans fc Dicke 



1990 



1961 



Dicke 



De Sabbata et al. 



1992 



Capozziello fc de RitisI 



Melnikov 



1994a 



Capozziello et al. 



1996 



Melnikov 1996 



19971: ICapozziello et al.lll998l: iDrinkwater et all Il998l: iFischbach fc Talmadgd Il998 



Barrow fc O'Tbolel l200l[ IPanielssonl I2001I: 



Krause fc FischbachI 2001: Murphv et al. 



2001 



Garcfa-Berro et al. 2007: Bertolami et al. 



Uzanl I2OO3I: Ide Sabbata et all 120041: 



2008 



Capozziellol 20051: Clifton et al. 



Hamber fc Williamsl I2OO5I: iBailev fc Kosteleckvl I2OO6I: iBrownstein fc Moffat 



Capozziello fc Francaviglia 



2OIOI : iKostelecky fc TassonI 



2OO5I: 



2006 



2008 



Uzanl 



Adelberger et all l2009l : lUzanI l2009l : iRami 



1976: Warburton fc Goodkind 1976 



Mikkelsen fc Newman 



1977 



2011 



and e xperimentally/ observationally ( Wagoner! 19 70l: Vintil 19721: Ulrichri974j:~Lond 



20111 ) 



Anderson et al 



Kisliklll983l: IChan fc Paik|ll984l: ICilliel 



1978 



Blinnikov 



1978 



Hut 



19811 : 



Chan et al. 



1982 



Burgess fc Cloutierl Il988l: iTalmadge et al.l Il988l: Krauss fc White! Il992l: llzmailov et al. 



1987 



1993 



Paik et al. 1994 



Bertolami fc Garcia-Bellido 1996 



Gillies 



19971: iGaztanaga et al 



2002 



LongI 



Gershtevn et al. 



2002 



20031: lAbraravan 



2009 



2004' 



2004 



Barrow 



Garcfa-Berro et al. I2OO7I : Ilorioll2007l: 



Unnikrishnan fc Gillied 2002alJbl: lAdelberger et al 



Lammerzah]||201ll: iPiedipalumbo et al. 



Milnd (119351, Il937t ). iDirad (119371 ) and IJordanI fll937l . ll939D . 



2011 



Uzan 



2003 



2OO5I: iBoucheijlgOgsl : jKruonogov fc Mernikovl booi : 



Bertolami fc SantosI l2009l: iNewman et all l2009l: [Li 



20111 ) since the early insights by 



In this paper we deal with possible smooth anisotropic spatial variations of G, i.e. we 
consider the case G = G{r) from a purely phenomenological point of view. We stress that 
in our analysis we do not rely upon any specific theoretical scheme encompassing such a 
spatial variability of G: the interested reader may want to consult the previously cited 
specialized literature. Quite generally, we express a putative anisotropic dependence of G 
on the spatial coordinates by parameterizing it with a gradient VG along a fixed direction 
in space ^ as 

G{r) ~ Go + VGo • r, (1) 



with 



VGo = I VGo I e 



(2) 



The subscript "0" in eq. ([T]) and eq. ([2]) refers to quantities evaluated at the origin of the 
spatial coordinates which, in our case, coincides with a generic body of mass M acting as 
source of the gravitational field. By means of eq. ([1]) we are assuming that the putative 
variations of G are rather smooth over the spatial extensions considered. A change of G 
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like that of eq. ([T]) is usually absent in the standard alternative metri c theories of gravity 
treated within the Parameterized Post-Newtonian (PPN) framework (jWilllll993l ) in which 
G may depend on the velocit y V of the frame in which the experiments are performed with 
respec t to a preferred frame ( 1Willlll97ll ): see, e.g., IVintil (119721 ) ■ iDamour fc Esposito-Faresd 
( I1994J ) and IWilll ( 1l993l )for analyses of the orbital consequences of such kind of anisotropies. 
A spatial anisotropy of G depending on the angle between the line of interaction of two 
gravitating bodi es and a reference direct i on wi th respect to distant stars was experimentally 
investigated bjfl iGershtevn et"aD fcool l2004l ) in a series of Earth-based laboratory 
i nvestig a tions which were subsequently critically analyzed by lUnnikrishnan fc Gillies 



(|2002aJ). lUnnikrishnan fc GilliesI ( l2002a( ) remarked that, in general, spatial anisotropies 



of G may occur depending on how nearby masses and their distribution can affect the 
gravitational interaction between two bodies, and also because of preferred frame effects. 
If the modification of the gr avitational interaction depends on the gravitational potential 
generated by other masses ( iBrans &: Dickd Il96ll : ISciamal Il953[ ) in a somewhat Machian 
fashio n, then the most distant ones d ominate and the spatial anisotropy is expected to be 
small ( lUnnikrishnan fc GilliesI l2002al ). It is, then, reasonable to expect that the Galaxy 
yields the most important contribution to the anisotropy by assuming ^ directed towards 
the Galactic Center (GC) (lUnnikrishnan fc GilliesI l2002al ). 



In the framework of our parameterisation of eq. ([T]) and eq. (|2]) , such a scenario may 
offer, in principle, interesting observational perspectives if suitable astronomical bodies are 
chosen. Indeed, if we insert eq. ([1]) in the usual expression of the Newtonian inverse-square 
law, a small modification of it occurs 



A 



M(VGn-r)r M\VGo\U-r]r 



Now, the ecliptic coordinates of the GC are (IReid fc Brunthalerll2004l ) 

Ago = 183.15 deg. 



(3) 



(4) 



/3. 



GC 



-5.61 deg. 



so that the angle between ^ and r for, say, a typical planet of the solar system is rather 
small. Thus, 

AG |VGo|r 



G 



Gn 



(5) 



Unnikrishnan fc GilliesI ( l2002al ) argued that the order of magnitude of the Galactic-induced 
anisotropy is 

(6) 



G 



^They claimed to have measured an anisotropy as large as 5.4 x 10 ^ with a torsion 
balance. 
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where c is the speed of hght in v acuum, Mpai ^ lO^^M fT) (IBattagha et al.ll2005[ ) is the mass 
of the Galaxy, and d = 8.28 kpc (IGillessen et al.l 120091 ) is the distance from the GC. In our 
picture, it would naively be equivalent to perturbing accelerations 



where the standard inverse- square Newtonian accelerations An are 



(7) 



4 X 10"^ m s"^ < An < 4 X 10" 



m s 



for the major bodies of the solar system. In principle, perturbing accelerations as large as 

A ~ 10"® - 10"^^ m s"2 (9) 



may have interesting observational consequences. We will analytically work out them in 
detail. Indeed, relying upon simple order-of magnitude evaluations may be misleading 
since important factors of the order of 0{e^), j = 1,2, ... or 0{e~^), j = 1,2, ... in 
the usually small eccen t ricities e of the bodies adopted as probes may be neglected. 
Unnikrishnan fc GilliesI f l2002a( ) performed a preliminary calculation concerning the Earth 

12 



Moon system. They started from a certain value 



G/G 



< 4 X 10"^^ yr"^ fiDickev et al 



19941 ) of the upper bound in the fracti onal time change of G obtained with the Lunar 
Laser ranging (LLR) technique. Then, lUnnikrishnan Sz GilliesI (j2002al ) stated that the 
same analysis can be useful as far as the spatial anisotropy of G is concerned as well. 
They noticed that an anisotropic spatial variation of G like that of eq. ([T]) should exhibit 
a harmonic signal with the same approximate monthly periodicity of the orbital lunar 
motion as the line joining the Earth and the Moon sweeps out different directions with 
respect to the GC direction. Finally, since the accuracy with which it is possible to measure 
a periodic signal may be of th e same order of, or better than that for a secular trend, 
Unnikrishnan fc Gillies (2002a) concluded by inferring 



AG 
~G' 



< 4 X 10 



-12 



(10) 



Repeating the same reasonings with the latest results from LLR (Williams et al. 20041 : 



Miiller fc Biskupekl 120071 : IWilhams et al.ll2009l : iHofmann et al.ll2010h would yield 



AG 
IT 



< (1 - 0.4) X 10^ 



-12 



'in 



in neat disagreement the results by iGershteyn et al.l (120021 ). iGershteyn et al.l ( l2002l ) 



pointed out that such a disagreement exists if it is assumed that the G anisotropy depends 
neither on the magnitude of the interacting masses nor on the distance between them; 



^It is so also because 1 year contains almost 12 lunar cycles. 
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Gershtevn et al.l (|2002[) remarked t hat the masses and the distances invo lved in the analys is 



bv lUnnikrishnan fc GilliesI (j2002al ) drastically differ from those used by iGershteyn et al. 
( I2OO2I ) in their experiment. 



We propose to obtain much more accurate bounds than that in eq. (fTOjl by calculating 
in detail all the orbital effects of eq. ([3]) which, actually, depends on the mutual distance 



far as other performed and/or proposed astronomical tests of G(r) are concerned ( 


Waeoner 


197nl: Vinti 19721: Ulrich 1974: Warburton & Goodkind 


19761; 


Vdikkelsen & NewmanI 1977: 


Anderson et al. 


1978 


Blinnikovl 


1978|; 


Hut| 


1981; 


Kislik 


1983; 


Bureess & Cloutierl 


19881: 


Talmadee et al. 


1988 


Abramvan 


2004; 


lorio 


2007; 


Bertolami & Santos 2009: Li 2009f). quite 



different explicit theoretical isotropic models and empirical approaches have been followed 
so far; in particular, an exponential Yukawa-type model and the third Kepler law have often 
been adopted. 

The paper is organized as follows. In Section [2] we analytically work out various 
orbital effects caused by a phenomenological dipole-type spatial variation of G like that of 
eq. (II]) averaged over one orbital revolution of a test particle. In Section |3] a comparison 
with the latest solar system planetary observations is made; we also consider the stellar 
system around the Galactic black hole. Section S] summarizes our findings and contain the 
conclusions. 



2. Calculation of the orbital effects 



The orbital motions of test particles are, in principle, a ffected by eq. (|3l) whose eff ects 



can be worked out with standard perturbative techniques (iBrouwer fc Clemencd 11961 
Bertotti et aPboOSh . By defining 



Grn 



(12) 



the radial component Ar of eq. ([3]), evaluated onto the unperturbed Keplerian ellipse, is 



A 



R 



tpoHo (1 + e cos /) |cos u (^i^^ cos i? + sin J? j + sin u sin / + cos / (^i,y cos Q — sin i? 



a(l-e2) 



(13) 

the transverse and out-of-plane components At-, A]^ vanish. In eq. f lT3|) a, I, f2,u = u+f, u, f 
are the osculating semi-major axis, the inclination of the orbital plane to the reference 



^Notice that [tpo] = while [GoM] = T'^ as usual. 
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{x, y} plane, the longitude of the ascending nod^, the argument of latitude, the argument 
of pericentral and the true anomal}|§, respectively, of the test particle. Notice that eq. ( !T3|l 
clearly shows that the putative anisotropic G effect depends, among other things, on the 
distance between the two bodies as well. By assuming ^ constant over one orbital revolution 
of the test particle, a straightforward first-order applicati on of the Gauss perturbative 
equations ( Brouwer fc Clemence 1961 : Bertotti et al. 2003 ) yields for a, e, the longitud^ of 
the pericenter w = f2 + u, and Ai the following non-vanishing long-term temporal rates of 
change 



da 

dt ~ 



de 
dt 



aen 



I — sin cos + iy sin i? j -|- cos u sin I + cos / {^y cos Q — S,x sin i? j | 



i/'OMo(l-e2) 



— sin u ( ^x cos i? -|- sin i? j -|- cos u 



sin / -|- cos / ( cos Q — sin i? 



dvj 
dt 



i/'OMo(l-e^) 



COS u ( COS Q + ^y sin i? ) -|- sin uj 



sin / -|- cos / I ^y cos Q — ^x sin i? 



dM 

dt 



'/'OMo(l-e^) 



3/2 



COS uj ( ^x COS Q + ^y sin i? ) -|- sin uj 



^z sin / + cos I y^y cos Q — ^x sin Q 

(14) 

The inclination J and the node Q are left unaffected. Notice that the quantity n entering 
eq. (JH]) is the unperturbed Keplerian mean motion, i.e. n = ^^/a^. The long-term rates 
of change of eq. (fT^ are exact in the sense that no a-priori assumptions concerning ^ and 
the orbital configuration of the test particle were made. All the formulas of eq. (fT^ are 
singular for e — )■ 0; however, it is just a n unphysical artifact which can be cured by adoptin g 
the well-known non-singular elements (IBrouwer fc Clemencdll96ll : iBroucke fc Cefolalll972l ) 



h 



esmttJ, 



ecosccj. 



(15) 



/ = w + M. 



■^It is an angle in the reference {xy} plane counted from the reference x direction to the 
line of the nodes, which is the intersection of the orbital plane with the reference {xy} plane. 

^It is an angle in the orbital plane reckoning the point of closest approach with respect 
to the line of the nodes. 

^It is a time-dependent angle in the orbital plane determining the instantaneous position 
of the test particle with respect to the pericentre. 

^It is a "dogleg" angle. 
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It is also important to remark that the vahdity of eq. f lMj) is not restricted to any specific 
reference frame, being, instead, quite generaL 

It is useful to work out the radial, transverse, and out-of-plane shifts over one orbital 
revolution of the positio n and ve l ocity vectors r and v as well. They can be analytically 
worked out according to ICasottd (119931 ). For the shifts Ai?, AT, AA^ of r we have 



Ai? 



— sin w ( cos i? + sin i? j + cos u 



sin / + cos / ( cos H — sin i? 



AT = ^""^0° ^ J ^ ^ / j-.Qg ]7 _|_ gi]2 j7 J -|- sin uj C,z sin / + cos / ( C,y cos i? — sin J? 



AN 



0. 



The shifts Avji, Avt, Ati^r of v are 



(16) 



Avr 



27rV)oa^n(l+e)^/^ 



Vl— e 



COS cu ( .^x- COS i? + sin i? ) + sin cu 



sin J + cos I ( .^y cos Q — sin 12 



Afy 



27rt/ioa'^n\/T^ 



sin UJ { ^x cos + sin i? j + cos cu ^z sin / + cos / ( cos Q ~ ^x sin i? 



Awtv 



0. 



(17) 

Concerning the singularities occurring for e — 0, also for eq. ( IT6|) and eq. (fT7j) the same 
considerations as for eq. f|T^ hold. 



Confrontation with the observations 



Here we compare our theoretical results of Section [2] with the latest observationally 
determined quantities pertaining the orbital motions of some of the major bodies of the solar 



(Pitjeva 


2007; 


Fienea et al. 


2011. 


2010; 


Pitjeva 


2010) 



bounds on ipQ. In principle, one should explicitly include eq. ^ in the dynamical force 
models fit to the observations and re-process the entire data set with such an ad — hoc 
modified theory by varying the parameters to be estimated, the data, etc. Otherwise, the 
putative signal may be partly or totally absorbed in the estimation of, say, the initial state 
vectors. However, this is beyond the scopes of our work. 

In computing the anomalous effects for different bodies, we refer their orbital 
configurations to a heliocentric frame with mean ecliptic and equinox at the epoch J2000. 
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In it the unit vector pointing to the GC is 

= -0.993, 



-0.054, 



-0.097. 



In Tabled] we quote reasonable evaluations for the secular variations of the semi- major axes 
a and the eccentricities a of the inner planets of the solar system for which the most accurate 
data are currently available. They were co mputed by div iding the formal, statistical 1 — cr 
errors in a, e of the EPM2006 ephemerides iPitjeval ( 120071 ) by the time interval At = 93 yr 
covered by the observations used for constructing them. The realistic uncertainties may 



Table 1: Formal uncertainties in the secular variations of the semi-major axes a and the 
eccentricities e of the inner planets of the solar system. They were obtained by dividing 
the formal errors in a and e in Table 3 of iPitjeval (120071 ) by the time interval At = 93 yr 
(1913-2006) of the data records used for the EPM2006 ephemerides (IPitieval lioOTh . The 

errors in e were computed as ae = sj 
up to one order of magnitude larger. 



\dh> 



0"?. The realistic uncertainties may be 



Planet 


(Ta (m yr ^) 


(^e (yr ^) 


Mercury 


3 X 10-3 


4.43 X 10-12 


Venus 


2 X 10-3 


1.8 X 10-13 


Earth 


1 X 10-3 


5 X 10-1^ 


Mars 


3 X 10-3 


5 X 10-1^ 



be up to one order of magnitude larger. Table [2] displays the latest determinations of the 
corrections Aw to the standard Newtonian/Einsteinian secular precessions of the perihelia 
of th e planets of the solar system recently obtained with different ephemerides (IFienga et al. 
201 ll . I2OIOI : IPitjeval |2010[ ). Since ipQ enters eq. (fT^ as a free, adjustable parameter, we 
select the values for it which make the putative anomalous effects of eq. (IHl) not larger 
than the empirically obtained bounds in Tableland Table [2l As a result, we are able to 
obtain reasonable guesses concerning AG/G at different heliocentric distances by posing 
AG ~ (T^p (r). The results are shown in Table [31 Even by re-scaling the bounds obtained 
from Tabled] by one order of magnitude, the anisotropy of G in the planetary regions of the 
solar system is very tightly constrained, being of the order of 10-^^ — 10-^''. 

It may be interesting to consider a completely different astronomical scenario, both 
from the point of view of its components and of the distance scales involved. In Table H] 
we quote the relevant physical and orbital parameters of the system constituted by the 
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Table 2: Estimated corrections Acij, in milliarcseconds per century (mas cty^^), to the stan- 
dard Newtonian-Einsteinian secular precessions of the longitudes of the peri helia "^^7 of th e 
eight planets of the solar system p lus Pluto determined with the EPM2008 (|Pitievall2010l ). 
the INPOP08 fiFienga et al.ll2010h . and the INPOPlOa flPienga et all l201lh ephemerides. 
Only the usual Newtonian-Einsteinian dynamics was modelled, so that, in principle, the 
corrections Aro may account for any other unmodelled/mismodelled dynamical effect. Con- 
cerning the values quoted in the third column from the left, they correspond to the smallest 
uncertainties reported by iFienga et al. fl2010f ). Note the small uncertainty in the correction 
to the precession o f the terrestrial perihelion, obtained by processing Jupiter VLBI data 
fiFienga et al.ll2010h . 



Planet 


fPitieva 2010) 


Attj fPienga et al. 2010) 


Aro fPienga et al. 2011) 


Mercury 


-4 ±5 


-10 ±30 


0.2 ±3 


Venus 


24 ±33 


-4 ±6 




Earth 


6±7 


0± 0.016 




Mars 


-7 ±7 


0±0.2 




Jupiter 


67 ±93 


142 ± 156 




Saturn 


-10± 15 


-10 ±8 


±2 


Uranus 


-3890 ± 3900 


± 20000 




Neptune 


-4440 ± 5400 


± 20000 




Pluto 


2840 ± 4510 







Table 3: Upper bounds on the anisotropic percent variation AG/G inferred from d, e, for 
the inner planets of the solar system. We posed AG < cr^;, (r) = o^^a (1 ± e^/2) for each 
planet, where was obtained by comparing the theoretical predictions of eq. f[T^ to the 
uncertainties listed in Table [T] and Table [H 



Planet 


AG 
G 


a 




AG 
G 


e 




AG 
G 




Mercury 


2 X 10" 


-16 


8 X 10" 


-15 


3 X 10- 


-13 


Venus 


9 X 10- 


-18 


1 X 


10- 


-18 


1 X 


10- 


-17 


Earth 


1 X 


10- 


-17 


2 X 10- 


-18 


3 X 10- 


-18 


Mars 


4 X 10- 


-16 


3 X 10- 


-16 


3 X 10- 


-15 



Supermassive Black Hole (SBH) hosted by the GC in Sgr A* and t he main sequence sta r 
S2 orbiting it in about 16 yr at a distance of approximately 1 kau (IGillessen et al.ll2009l ). 
In this case, the angular elements refer to a coordinate system whose reference z axis is 
directed along the line of sight: the reference {x^y\ plane coincides with the plane of the 
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sk}|§, with the x axis pointing towards the Celestial North Pole. The results of eq. 



Table 4: Relevant physical and orbital parameters of the SBH- S2 system in Sgr A* . The 
Keplerian orbital elements of S2 were retrieved from Table 1 of (IGillessen et al.ll2009[ ). The 
figure f or the gravitational parameter comes from a multi-star fit yielding fiQ = 4.30 x 



10^/^0 (IGillessen et al.ll2009l ): it yields a Schwarzschild radius as large as Vg = 0.084 au. The 



quoted value in m for the semi-major axis of S2 was obtained by multiplying its angular 
va lue a = 0.1246 arc sec ( IGillessen et al.l 120091 ) by the distance of the SBH d = 8.28 kpc 



(IGillessen et al.ll2009l ): it corresponds to 1031.69 au, so that the orbital period of S2 is 
Pb = 15.98 yr. The uncertainty in the secula r precession of S2 ca n naively be obtained 
by dividing the error in u, quoted in Table 1 of ( IGillessen et al.l 120091 ). by the time interval 
AT covered by the observations used which is almost equal to Pb- 



fiQ (m^ s ^) a (m) 



/ (deg) f2 (deg) u (deg) cr^j (arcsec yr ^] 



5.70 X 10^6 1.54 X 10^^ 0.8831 134.87 226.53 64.98 



182 



are applicable to S2 as well since its mass is about five orders of magnitude smaller than 
that of the SBH: clearly, in this case it is ^ = —z. Concerning the use of eq. ( !T^ in 
the Sgr A* scenario, one may wonder why an essentially Newtonian scheme is adopted 
instead of a general re lativistic one. In principl e, one could assume as reference path a fully 
post-Newtonian on^ ( Calura et al. 1997 , 1998), and work out the effects of a given small 
extra- acceleration like eg. ([3l) w ith respect to it according to the perturbative scheme set 



up by ICalura et al.l (Il997l Il998l ). which is a general relativistic generalizatio n of another 



stand ard perturbative approach based on the planetary Lagrange equations ( iBertotti et al. 



20031 ). Actually, it is, in practice, useless since the only addition with respect to the orbital 
effects like, e.g., the precession of the pericenter, resulting from the standard scenario would 
consist of further, small mixed GTR-perturbation orbital effects, completely irrelevant in 
strengthening the bounds inferred. Viewed from a different point of view, the ratio of the 
average distance rs2 of S2 from the SBH to its Schwarzschild radius Vg is, after all, as large 
as 1.7 X 10^ 

By using the expression for the putative precession of the stellar pericenter in eq. (IT4|) . 
it can be compared to the present-day uncertainty in observationally determining its secular 
rate in Table |H In this case, the constraints on the G anisotropy are very weak. Indeed, we 



'^It is a plane tangent to the Celestial Sphere at the point where the object of interest is 
located. 

^In doing so, it would be implicitly assumed that the effects due to a putative G anisotropy 
are smaller than the post-Newtonian ones as well. 
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have just 
corresponding to 



^"(7^ > 0.45 pc, 
AG 



< 1.5 X 10" 



(19) 



(20) 



over about 1 kau. Note that the bound of eq. (|2U]) is tighter by one order of magnitude 
than that could be inferred by simply posing 



— <^ = 1.2xl0-\ 
G /io 



(21) 



from = 0.50 x 10 Vo flGillessen et al.ll2009h . 



4. Summary and conclusions 

We looked at phenomenological anisotropic spatial variations AG/G of the Newtonian 
gravitational coupling parameter G, in the form AG = VG • r, and analytically worked 
out the impact that they may have on the trajectory of a test particle orbiting a central 
body of mass M. More specifically, we focussed on the cumulative orbital changes obtained 
perturbatively by averaging over one period of revolution of the test particle the effects 
due to AG/G on its path. As a result, the osculating semi-major axis a, the eccentricity 
e, the pericenter zu and the mean anomaly Ai of the orbiter experience no n- vanishing 
long-term changes which depend on the overall orbital geometry of the test particle and 
on the direction ^ of the putative gradient VG. We analytically worked out the long-term 
variations per orbit AR, AT, AN and Avji, Avt, Av^ of the position and velocity vectors 
r and v of the test particle as well. We found that both the radial and the transverse 
components of r and v are affected by long-term changes per orbit, while the out-of-plane 
ones are left unaffected. We kept ^ fixed during the integrations: moreover, no a-priori 
simplifying assumptions on e and I were assumed, so that our results are exact in this 
respect. 

Then, we compared our theoretical predictions to the most recently determined 
observational quantities for some of the major bodies of the solar system. By assuming that 
the dominant contribution to the hypothetical anisotropy of G is due to the Galaxy, we 
took ^ directed towards the Galactic Center, which has a small inclination with respect to 
the ecliptic. By using the heliocentric orbits of the inner planets we were able to constrain 
AG/G to a level of about 10~^^ over ~ 1 au, several orders of magnitude better than 
in previous analyses based on Lunar Laser Ranging only. We looked also at the star S2 
orbiting the Galactic black hole at a distance of about 1 kau along a highly elliptical ellipse, 
but, in this case, we got just AG/G < 10"^. 
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